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Abstract 

We examine the problem of determining which nodes are neighbors of a given one in a 
wireless network. We consider an unsupervised network operating on a frequency-flat Gaus- 
sian channel, where K + 1 nodes associate their identities to nonorthogonal signatures, trans- 
mitted at random times, synchronously, and independently. A number of neighbor-discovery 
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algorithms, based on different optimization criteria, are introduced and analyzed. Numerical 
results show how reduced-complexity algorithms can achieve a satisfactory performance. 

1 Introduction 

Of late, wireless networks, and in particular sensor networks, have been the object of a good deal 
of interest, also spurred by the manifold applications they can be associated with (see, for exam- 
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pie, their applications to classification and tracking [1] and to monitoring [2]). A characteristic 
requirement of several wireless networks, which enables them to adapt themselves to a changing 
environment, is that they be "self-configuring," i.e., that a large number of wireless nodes orga- 
nize themselves to perform the tasks required by the application they have been deployed for: 
examples of self-configuration include construction of routing paths, clustering, and formation 
of minimum- weight trees. In this paper, we consider an aspect of self-configuration in wireless 
networks referred to as neighbor discovery (ND). Neighbor discovery is the determination of all 
nodes in the network a given node may directly communicate with. Knowledge of neighbors is 
essential for all routing protocols, medium-access control protocols, and several other topology- 
control algorithms. Ideally, nodes should discover their neighbors as quickly as possible, which 
will allow nodes to save energy in their discovery phase. Also, rapid discovery allows for other 
protocols (such as routing protocols) to quickly start their execution. In addition, ND may also be 
the solution for "partner selection" in cooperative wireless networks. In fact, cooperation among 
users may carry advantages only if the partners are chosen in a proper way: for example, "decode- 
and-forward" (DAF) protocols may suffer from cooperation with weak users, thus failing in the 
goal of increasing the diversity order [3]. 

Recently, a number studies on ND algorithms have appeared (see, e.g., [4, 5] and the refer- 
ences therein). Most of these approach ND at a protocol level, defining node A to be a neighbor 
of node B if A can exceed B's signal to noise-ratio requirement: as a consequence, A is inserted 
in the neighbor list of B based solely upon successful reception, at node B, of a packet sent by 
node A. Moreover, the Internet Engineering Task Force proposes to perform Neighbor Discovery 
"at IP Layer" [6]. The corresponding protocol assumes a broadcast capability at physical layer, 
and a MAC which handles contention. Now, ND algorithms for wireless networks may not be 
contention-based when energy constraints are tight: retransmission in the case of a collision costs 
energy, which might be a resource at a premium. In this context, we consider a transmission 
scheme which avoids collisions at modulation level and is based on simultaneous transmission of 
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signatures. In principle, if the nodes' waveforms were orthogonal, no collision would occur. In 
practice, these waveforms exhibit a small correlation, which causes an interference, red In a "stan- 
dard" network approach, signature transmission fails whenever there is interference, because this 
causes a "collision" which may occur very often during a ND session if there are many neighbors. 
Using multiuser detection, interference does not cause collisions (these are in a sense "automati- 
cally resolved") and can be controlled by multiuser-detection algorithms. 

ND can be performed in a supervised or unsupervised manner. In supervised methods, there is 
a central controller (e.g., a leader node) which processes the signal received from all nodes, deter- 
mines the network configuration, and communicates to all nodes their neighbor lists. Supervised 
ND algorithms are expected to cost a large amount of energy, and hence they should be discarded 
for energy-limited networks, like sensor networks. Unsupervised ND algorithms have no central 
controller, and each node discovers its own neighbors. Another important issue in ND problems 
is the timing aspect. In [7], the frame-synchronous assumption is justified by the presence in each 
node of Global Positioning System (GPS) devices. In [5], asynchronous algorithms are addressed, 
assuming that nodes can synchronize at bit level (which is the assumption we make in the follow- 
ing). 

The goal of this work is to provide the foundations of signal processing for ND in wireless 

networks. We consider an unsupervised wireless network, where ND is operated independently of 

the regular exchange of packets in a frequency-flat Gaussian multiple-access channel, shared by 

K+l nodes which transmit, synchronously and independently, a set of known signatures according 

to the scheme advocated in [5]. Each node is identified by its own unique signature, and every node 

keeps a list of all the signatures of the network[] A node is called a neighbor of the reference node 

'We do assume that ND can occur in a separate channel of a mobile system. In a sensor network, where energy 
is at a premium, it makes sense to perform ND as soon as nodes are deployed. Other ND algorithms can be based 
on higher layer protocols such as IP, but these might seriously waste energy. Even if TCP/IP is in use (in a mobile 
network, say), then allocating a small fraction of overall system resources at the physical layer to ND might be better 
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if its amplitude, received by the latter, exceeds a preassigned activity threshold, say tau Moreover, 
nodes cannot transmit and receive simultaneously on the same channel^] and the maximum number 
of active nodes is fixed and finite. We clarify that a neighbor relation between two nodes need not 
be bidirectional, since each node discovers those nodes it can receive from. 

The organization of this paper is the following. In Section[2]we provide a model for the physical 
aspects of the networks, and we formulate our problem. ND algorithms are introduced in Section[3l 
and analyzed in SectionSl Section |5] shows some numerical results, while Section[6]concludes the 
paper. 



2 Signal model and problem formulation 

Our scenario is based on the transmission scheme illustrated in Fig. [T] which corresponds to node 
searching its own neighbors among four other nodes. Q In every time interval ("slot"), each node 
i, i = 0, 1, . . . , K, transmits its own signature, independently of the other nodes, with probability 
Si, while otherwise (and hence with probability 1 — £j) it senses the channel. This probability is 
actually designed as a part of the algorithm, and it influences the ND algorithm performance, as 
we shall examine in our analysis. 

The ND algorithm runs in a finite period, called a discovery session, whose duration is denoted 
To- During Tp, every active node transmits a number of signals containing one or more copies of 

than running ND periodically as an application program. 

2 Note that this definition can be generalized: for example, one may define a neighbor as one whose power-to- 
interference plus noise ratio exceeds a given threshold. In this paper we stick to a more restrictive definition, which 
allows simpler algorithms. In general, any information gleaned through ND will help scheduling, or medium-access 

control, or routing algorithms to be more efficient, because they will be armed with neighbor information. 

3 For simplicity, we disregard the more general case of nodes that can be in an idle state, i.e., they are neither 

receiving nor transmitting. 

4 We consider node to be the reference node. Since all nodes are at the same hierarchical level, the same analysis 
applies to any node. 
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its signature. Each signal has duration T = Tjj/N, with N the number of slots in the discovery 
session. The network is assumed to be unsupervised, which implies that all nodes are independent 
and at the same hierarchical level: as a consequence, the ND algorithm is run in parallel by all 
nodes. Under the assumptions made in Section[Q the baseband representation of the signal received 
by node in the time interval \(n — 1)T, nT\ n G {1,2,..., N}, is 

]Cf=i ^k, n a k Sk (t- (n- 1)T) + z(t) if tp , n = 

y[t) = < (i) 

[ ifVv = ! 

where denotes the channel gain, i.e., the complex amplitude of the signal received from node k 
and assumed to be constant during all the discovery session, s k (-) is the A;th node signature, vjjk, n 
is a random variable taking value 1 if node k is transmitting at time n, and value otherwise (so 
that F(ijjk,n = 1) = and z(t) is additive white complex Gaussian noise having spectral density 
2N . We assume a>k to be modeled by a complex circularly symmetric Gaussian random variable 
with variance 2a\. The signatures can be expressed as 

L 

Sk(t) = Sl M f - V - l ) T c)/^L (2) 
i=i 

where G {—1, +1} is the Zth chip of the A;th signature, L is the processing gain, T c = T/L 
is the chip duration, and </>(•) is the (unit-energy) chip waveform|f| The slots devoted to channel 
sensing need not be adjacent: however, due to our flat-fading assumption, we may assume, without 
any loss of generality, a sensing phase of 

N 

M = Y / (l-4 0)n )=N-u (3) 

71=1 

consecutive slots with intermittent other-users activity, with u the number of slots where node "0" 
is transmitting. Notice that M is random (N is assumed fixed and node has its own activity 
factor 6o), but the value it takes is known to node 0. Hence, in all subsequent derivations we refer 
5 The signatures are assumed to have unit energy. 
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to a given value of M . Of course, we may adopt the silent phases of node as a time scale, 
recasting (OQ), with a slight notational abuse, in the form£ 

K 

y{t) = J2^k, P a k s k (t-(p-l)T)+z{t), < t < M T, p=l,2,...,M (4) 

k=l 

Our problem is now reduced to determining the indexes k such that {|a/ c |}^L 1 exceed an "activity 
threshold" ta, based on model dU). 

Since z(t) is white Gaussian noise, the components of y(t) orthogonal to the subspace spanned 
by the signatures are irrelevant to our detection problem [8]. As a consequence, we might in 
principle adopt the signatures themselves, and their delayed versions, as an expansion basis for 
such a subspace. Alternatively, we may use the L— dimensional orthonormal basis 

L-l 

U{0(t-£T c -(p-l)T)} (5) 

e=o 

to expand the signal in the interval [(p — 1)T, pT) . The two approaches are obviously equivalent, 
but the latter is mandatory in situations where the discovering node has no prior information as to 
the signatures of other users: although we do not deal blind ND in this paper, we choose this one 
due to its inherent flexibility. 
Defining the scalar products 

rpT 
'(p-l)T 

with * denoting conjugation, we obtain a vector representation y p = [yi tP , y2, P , ■ ■ ■ , 2/L.p] T [] of the 
signal received in [(p — 1)T, pT) : 

K 

y P = ^2 ^ k >p akSk + z p = S* p a + Zp (7) 



(■pT 

Yi, P = / y(tW (t-(i- l)Tc - (p - 1)T) dt (6) 

J(v-1)T 



k=l 



6 Notice that the index n refers to consecutive time slots, while p refers to the time scale defined by the silent phase 

of node "0". 

7 The symbol T denotes transposition operation 
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where s fc = ^[sx, k , s 2 ,k, • • • , SL,k\ T , S = [s 1; s 2 , . . . , s K ], * p = diag(^i, p , ip 2 , p , • • • , ^k, p ), ol = 
[oil, oi 2 ,..., ot K ] T , z p = [z hp , z 2 , p , z L J T , and 

/ z(t)<j>* {t-(i- 1)T C - (p - 1)T) dt (8) 

J(v-1)T 



A 

Z i,'P 



'(p-l)T 

The ND problem now consists of assessing, after observing the set of M vectors yi:M — {yi> • • • ,yM }> 
which ones, among |aii|, . . . , \olk\, exceed the "activity threshold" ta- 

3 ND algorithms 

A sensible criterion for the selection of a ND algorithm consists of minimizing the probability of 
choosing, among the K network nodes under scrutiny, an erroneous set of neighbors of node 0. 
Since there are 2 K such sets, each corresponding to one hypothesis H, this error probability is 
minimized by the maximum a posteriori (MAP) decision rule: 

H = arg max P(H)p(y 1:Mo | H) (9) 

H 

where P(H) is the a priori probability of hypothesis H, and p(y~i : M \H) is the probability density 
of the observations given H. As an example, if K = 2, the 4 hypotheses are shown in Table [T] 

Hi H 2 H3 H4 

\ a i\ < r A I ckx I > ta, \ol\\ < ta |«i| > ta 
\a 2 \ < ta \a 2 \ < ta \a 2 \ > ta \a 2 \ > ta 



Table 1 : Hypotheses on the set of neighbors of node in a network with 3 nodes. 

Now, p(yi-.M I H) depends on the actual pattern of transmit/receive intervals of each node, 
denoted Vf^Mo- Since this is unknown under our assumption that the transmission of signatures is 
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not coordinated, it should be obtained from the marginalization 

1:M I H, *i:M ) 

*1:M 

which has a complexity that grows exponentially with KM . 

To overcome this complexity obstacle, the decision on the neighbor set works as follows. We 
first obtain estimates of the instantaneous powers of all nodes, next we decide that a node is a 
neighbor by comparing each of them with a threshold, i.e., 

H 1 

N 2 < rt (io) 

Ho 

where 

Hi'. The received instantaneous power exceeds r\. 
H : The received instantaneous power is below t\. 

(i) 

The performance of this test can be expressed through its probability Pp' of a. false-alarm and 
its probability P$ of a miss, defined as: 

P{ |«j| 2 > T? | \oii\ < T A \ 

(ID 

Hi > T A \ 

m 



P|? = P{H 2 <7f |N>T.4} 



These are related to the overall error probability througl 

P«(e) = PjpF{\ai\ < r A } + P^Pfel > r A } (12) 

Now, the maximum-likelihood (ML) estimators of the instantaneous powers can be obtained 
by jointly estimating a and the matrix sequence ^i-m - Straightforward calculations show that 



8 In what follows, the superscripts will be skipped whenever no confusion is induced by this notational simplifica- 
tion. 
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the ML estimates of ol and &i : m result from the solution of the 2 KM ° linear systems — each 
corresponding to an outcome &i;M ,i °f me matrix sequence *1?i : m '- 

'M \ /M \ 

vP=l / \p=l / 

with ^ denoting Hermitian operation. Computing 

Mo 

olml = arg min /A\y p - S&pjCti || 2 (14) 

i=l,...,2 KM f 

p=i 

with the solution corresponding to V&^mo,^ and recalling that ML estimates commute under 
nonlinear transformations, test (flOl) can be implemented by using |«j| 2 = |a?ML,i| 2 > 

Even with this receiver, implementation complexity would be unrealistic, and hence a further 
simplification is called for. Instead of dealing with the receive/transmit pattern related to the whole 
discovery session, we rather obtain estimates based on a single T-interval observation, which are 
then combined according to a suitable integration strategy. 

3.1 Suboptimum ND algorithms 

Consider again model ©. The ML estimate of ty p ct, based upon the observation y p available in 
slot p, is 

%Ta = (S^S)- 1 S^y p = S + y p (15) 

where S + denotes the pseudo-inverse of the tall matrix S. 
A closer look at this solution reveals that, since 

S + y p = S + SV p (x + S + z p = + w p (16) 

with E[iu p i/;t] = 2N (S^S)~ 1 , the interference from the other users is completely eliminated, 
at the price of some noise enhancement, reflecting the increase of the variance of its ith compo- 
nent by the factor {(S , '<S , )^~ i 1 }^ 1 . It is interesting to notice that this estimate is noise-limited, but 
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not interference-limited, implying that any receiver based on (1151) is asymptotically efficient [10]; 
likewise, near-far resistance is granted [10]. 

Since there are M sensing phases, the M estimates resulting from repeated application of 
(fT5l) should be combined to yield the final test statistic. Borrowing techniques from radar detection 
theory, reasonable combination criteria are coherent integration (CI), wherein an estimate of the 
instantaneous power is obtained as 



\a 



*\ ci 



, M 



Mi 



P =i 



(17) 



and incoherent integration (II) 



, Mo 

^ 2 n = wY,\( s+ yp) 



(18) 



p=l 



Notice that 



E 
E 



\a 



*l ci 



a 



i\ ii 



M , |a 4 | 2 
M , \ ai \ 2 



Si + 



1 - e» 
M 



+ 



Mo 



(X; 



2 e, + 2N (&S)Z} 



(19) 
(20) 



implying that both |«i| 2 w and |aj| 2 c/ can be interpreted as biased estimators of the instantaneous 
power received in each slot from node i: biases can however be absorbed in the detection thresholds 
T;, while what matters here is that they are both consistent in the mean square sense, a property 
that will be exploited later on. Inserting (fT7l) and ([18]) into (flOl) . and skipping factors that can be 
absorbed in the detection thresholds, we obtain the coherent detector (CD) 

2 



and the Incoherent Detector (ID): 



E^\(S + y p ) 



I 2 > T- 2 



| 2 <T- 2 



node z is a neighbor 
node i is not a neighbor 

node z is a neighbor 
node z is not a neighbor 



(21) 



(22) 
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Notice how the CD can also be interpreted in a different way. Indeed, it may be obtained by 
first pre-processing the observations so as to form the cumulative sum: 

M M , K \ K 

y = ^2 = ( ^ fc .p afcSfc + z p ) = ^fcttfcSfc + z = SVa + z (23) 

p =i p =i V fc=1 / fc=1 

where 

M M 

and V = diag(Vi, . . . , i^), then multiplying the new observation by S 4 " and finally extracting the 
ith component to form the test statistic (|2T|) . Rewriting equation (|23T) in the form: 

y= z^OjS^ +y^ j i> k a k s k + jz^ (25) 

useful signal v fc ^ noise 
interference 

with z ~ W c (0, 2A^oMoIl), where I/, is the L x L identity matrix, the CD is easily seen to be a 
member of the family of linear ND tests (LNDT), wherein a decision on the proximity of user i is 
made based on the rule: 



|cjy| 2 < rl (26) 

Thus, the CD (I2TT) can be also interpreted as the zero-forcing (ZF) member of the family (|26l ), 
obtained as the unique solution to the constrained minimization problem: 

c L zf = argmin Ci E[|c] Y%=i a k u k s k \ 2 ] 

\A') 

C i.ZF S i = $ 



with (3 7^ 0, which yield 



c lyZ F = (II - SiSf) s, = (28) 



where Si is the L x (K — 1) matrix obtained skipping the i— th column from S and 7^ denotes 
the projector onto the orthogonal complement of the column span of Sj. For future reference we 



9 Notice from d28l l that the parameter (3 2 has been set to ^grgy^i 
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remind here that [10] 

|cU F s fe | 2 = I (29) 
[ |s^ Sl | 2 =|| Si,!. || 4 if A; = i 

where t ± denotes the projection of s, on the above orthogonal complement: needless to say, since 

II s^! || 2 = 1/[(S* S)^l], the noise power is enhanced by a factor (S*S)~l. 

The vector Cj can be designed according to a number of different criteria. For example, in [12] 

an LNDT based on conventional matched filtering (MF), i.e., assuming 

Ci,MF = Si (30) 

has been proposed and analyzed for ND [12]. MF is indeed simple, but it results into interference- 
limited performance, as we shall prove soon, nor does it retain the near-far resistance property 
granted by ML-based detectors. 

A possible alternative to the ZF criterion is offered by the minimum-mean-output-energy (MMOE) 
strategy, first introduced in [11], wherein the vector Cj is obtained as the unique solution to the fol- 
lowing constrained minimization problem: 



Ci,MMOE = arg min Ci E 

C l,MMOE S i = 1 



(3D 



namely: 



M yy s * 

Ci,MMOE = IT I - — (32) 

where M yy = J2k=i ^k-^Wkl^l + 2N M I L . Due to the invariance of the decision rule to any 
positive scaling of the test statistic, an equivalent detector relies upon setting 

Ci,MMOE = M yySi (33) 



It might be worth recalling here that, since 



lim M'lsi oc ViS, (34) 



N -*0 yy 
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MMOE is itself asymptotically efficient. Likewise, it retains the near-far resistance property since 
the projection direction c^mmoe tends to become orthogonal to those signatures whose ampli- 
tudes become increasingly large [11]. The advantage of ([33]) over ZF is that it easily lends itself to 
adaptive implementations in situations where the signatures of the active users are unknown. Even 
though we do not deal with adaptive ND in this paper, we anticipate that a number of reduced com- 
plexity algorithms, ranging from the O(L) -complex Least Mean Squares to the 0(L 2 ) -complex 
Recursive Least Squares, can be easily applied for adaptive MMOE implementation. 



4 Analysis 



From now on we assume that the node to be detected is node "1". Consider first the ID. The 
conditional false-alarm and miss probabilities in assessing the proximity of node 1 can be written 

as: 

P M = P(xi < ^i 2 |l«i| > r A ,*l:M ) (35) 
P F = P(xi > 7?||Ql| < T A ,*i : Mo) (36) 

with xi — 2~2t=i l(S + y P )i| 2 - Given |ai| and Vf^Mo* Xi ls noncentral chi-square distributed with 
2M degrees of freedom and parameters ^i|ai| 2 and a\ 1 = (S^S)^}^, implying 

P(Xi > t\ I M, Hf 1:Mo ) = Q Mo ( y^M^ ^\ (37) 

V °n,i cr n> ij 

where Qm {'i ') i s me Marcum function of order M . Using the series expansion of modified 
Bessel functions 

w=E jg: t+1) ™ 
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we obtain 

k 



1 > T l I Kl> *l:Mo) = e ^ Er=0 2*ife!r(A/o+fe) F + ^ 2^ J 

*i:M ) = i+fer 2Xo (irSr) « ( M o + k; Q (* + 1; ^(1 + 



i > if, I ai| > T4 

(39) 

where 



r(fc; x) ± I t^e-'dt , Q(k; x) = (40) 

are the upper incomplete Gamma function and its regularized version, respectively, while pi is the 
signal-to-noise ratio after decorrelation, i.e.: 

^ = ^= .^tVi (41) 



<i iV (^5)r! 



We thus obtain the conditional measure: 



£(1^) o(M + *;^-)g(fc + l;g(l + ^) (42) 



1 + ^ fc =0 



which should be averaged over z/ x to yield the conditional probability of a miss given M . Similar 
developments hold for P F , yielding 

k 



P(Xi > 7f,|ai| < r A |*i :Mo ) = — Yl irZ^A Q(M +k;-^-)p(k 



(43) 

where 



7 (*; s) ^ / t^Vdt , P(*; x) = (44) 



o 



are the lower incomplete Gamma function and its regularized version, respectively. Finally, from 
(l43l) we easily obtain: 



P F = , cschf41rf " lPl 1 g(M + *;4-V(* 

2(l + ^pi) V4af;^Vl + W V 2<J V 



+ l;-4(l + z/i Pl 



(45) 



Analysis 



14 



DA, EB, ML: Neighbor Discovery in Wireless Networks 



February 5, 2008 



Consider now the test family (1261) . Notice that, since 

K 

9i = cjy = ^icjsiai + z/ fc cjs k a fc + c[z (46) 

useful signal k—2 ^ 
interference+noise 

\gi | 2 is conditionally chi-square with two degrees of freedom, given a\, {ui}f =l and M , with 
non-centrality parameter |z/ic{siai| 2 = |ai| 2 i/fc{siSi^Ci and scale parameter 

K K 

S 2 (ci) ^ l^ ct iSk| 2 ^ + M iV || Cl || 2 = ^a^clWci + MoiVolldf (47) 

k=2 k=2 

whereby, reproducing the same steps leading to (l42j) and (1431) . we obtain: 

A 

M 

i + v? 

k=0 



(48) 



(49) 



where p eq represents the signal-to-interference-plus-noise ratio (SINR) at the output of the linear 
filter, i.e.: 

Relationships d48l) and (|49l) are quite reminiscent of (1421) and (1431 ), respectively, one major differ- 
ence being the dependency of the performance on v\p eq , rather than V]_p\. Of course, the quadratic 
factor in v\ stems from the fact that linear detectors operate on a coherent combination of the ob- 
servations, while ID combines the slot-by-slot estimates incoherently. Notice, however, that the 
above relationships represent conditional measures, given M (i.e., given v ) and {z/j}fl 1 . If the 
discovery session is long enough, so that the matrix sequence vE^Mn ma y exhibit its typical behav- 
ior, namely, if N(l — e ) ^> 1, then the v k 's tend in probability to M £fc, whereby the unconditional 
performances may be obtained by averaging the corresponding conditional measures on the typical 
set of values of {vu}k=i an d M only, implying: 
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-M ~N(l-e ); 

- v k ~ M e k = Ne k (l - e ). 

In this limiting situation, it is interesting to notice the relationship between the "cumulated" SNR's 
for ID and CD (i.e., the ZF of ((281) ). i.e. (see also (f29b and subsequent comments): 

2 (J \ V \ II S l X II 4 V \°\ 

Thus, in terms of cumulated signal-to-noise ratio and for large N, ID seems to be preferable to 
CD, even though a global superiority cannot be claimed due to the different forms assumed by the 
respective false-alarm and miss probabilities. 

So far no criterion has been given to select the decision threshold T\. Notice, however, that the 
consistency of the estimates (fT71) and (fT8l) allows devising the asymptotically optimum thresholds 
(those achieving minimum error probability for large N) from (fl9l) and (|20l) in the form: 

rl CD = N(l-£ )[eni[N(l-e )ei + (l-e 1 )]+2N (S^S)i i W (52) 
r 2 l ID = N(l-e ) [>rle 1 +2N (&S)£] (53) 

For short discovery sessions, and under known activity factors of nodes to be discovered, optimum 
detection thresholds can be obtained by evaluating numerically the unconditional error probability, 
and then determining the points where it has a minimum. 

5 Results 

We consider here a fully loaded network with K + 1 = 7, each node being assigned a length-7 
m-sequence. As in the previous section, we assume that node "0" has to decide on the proximity of 
node 1. Figure |2] assumes SNRx = ol/N = dB, N = 100, a power-controlled scenario wherein 
all nodes are received with the same average power, uniform activity factor (e k = s = 0.5), 
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and an activity threshold equal to the median of the fading amplitude distribution, i.e., such that 
P(|cki| > ta) = 0.5. The figure represents the pair P M , P F for the various receivers examined 
so far. Interestingly, "conventional" MF suffers from the presence of the other nodes even in this 
rather benign situation, while MMOE, ZF (which coincides with CD rule) take advantage of their 
asymptotic efficiency: from now on, conventional MF will not be considered any longer. The 
interested reader is deferred to [12] for a fairly thorough performance assessment. 

The reliability of the asymptotic approximation for long discovery sessions can be assessed 
through figures [3]-[5]for the CD, and through figures [6]-[8]for the ID, which refer to the same scenario 
as in Fig. [2} The curves of these figures represent: 

- The unconditional false alarm and miss probabilities obtained by simulation. 

- The same pair obtained by a semi-analytical method, i.e., by estimating the averages of their 

conditional counterparts. 

- The asymptotic approximation. 

From the plots, it is evident that the asymptotic approximation tends to overestimate the perfor- 
mances in the interesting region of low error probabilities, while coming closer and closer to the 
true performance as N increases: notice that the approximation is extremely tight for N = 500, a 
realistic value indeed in real applications, which, for eo = 0.5, corresponds to M ~ 250. How- 
ever, it should be kept in mind that, for larger activity factors of the discovering node, the minimum 
value of N for the asymptotic behavior to be reached inevitably increases. 

The validity of the approximations (|52l) and (1531) . yielding the asymptotically optimal thresh- 
olds for CD and ID, respectively, can be verified through Figs. |9] and [RE showing P(e) versus 
T\ for CD and ID, respectively, for some values of SNR l5 e=0.5, iV = 500 and ta such that 
P(|ai| > ta) = 0.5. These plots, obtained numerically, show that the error probability admits a 
unique minimum; moreover, the optimal values of Ti are surprisingly close to those resulting from 
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the asymptotic approximations in (l52l) and (1531) . It might be marginally worth noticing that such 
an optimum value is practically independent of SNR4 for CD, while being strongly tied to SNR4 
for ID. 

Fig-QTlis aimed at comparing CD and ID. It represents the error probability versus the signal- 
to-noise ratio SNRx using the optimal thresholds for both receivers, and assuming again e = 0.5, 
iV = 500, and t a as before. It is interesting to notice that CD outperforms ID for small signal-to- 
noise ratios, while ID is preferable for medium-to-large values of SNR4. 

6 Conclusions 

We have examined the problem of discovery which nodes are neighbors in a wireless network 
operating over a fading channel. The optimum Bayesian decision rule has been derived, showing 
that its complexity is practically prohibitive. Two suboptimum neighbor-discovery algorithms have 
been introduced, based on standard techniques of coherent and incoherent integration. We show 
how coherent integration may be viewed as a particular case of a family of algorithm akin to Linear 
Neighbor Discovery Tests (LNDT). Theoretical analysis allows one to understand the design of a 
system employing such algorithms according to constraints on error rate, signal-to-noise ratio and 
discovery session duration. Finally, algorithm optimization was considered, and formulas were 
derived for asymptotical optimum threshold. 
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Figure 1 : A scheme for synchronous neighbor discovery. 




Figure 2: Performance of various ND algorithms under perfect power control, 2 of = 2N = 1 
(SNR 1= dB), N = 100, fully-loaded network. 
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Figure 3: Performance of the CD for 2o\ = 2N = 1 (SNR^O dB), N = 100 
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Figure 4: Performance of the CD for 2aj = 2N = 1 (SNRi=0 dB), N = 300. 
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Figure 6: Performance of the ID for 2o\ = 2N = 1 (SNR^O dB), N = 100. 
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Figure 7: Performance of the ID for 2crf = 2N = 1 (SNR^O dB), N = 300 
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Figure 8: Performance of the ID for 2o\ = 2N = 1 (SNR^O dB), N = 500. 
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